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Abstract. Solar wind plasma parameters and the Interplanetary Magnetic Field (IMF) 
observed by the WIND spacecraft upstream of the bow shock were used as input to 
magnetohydrodynamic (MHD) simulations of two substorm events. The power deposited 
into the ionosphere due to electron precipitation was calculated both from VIS observations 
and from the simulations. 


1. Introduction 

Global MHD simulations have been used for more 
than 20 years to help understand the interaction of 
the solar wind with the Earth’s magnetosphere 
[e.g. Leboeuf et al, 1978; Lyon et al, 1980]. 
Initially MHD models used idealized solar wind 
conditions (e.g. with IMF constant and 
northward or southward) to identify the principal 
features of solar wind-magnetosphere coupling 
[e.g. Ogino and Walker, 1984]. Advances in 
computers and improved niunerical techniques 
have allowed the construction of high-resolution 
three-dimensional models to study this interaction. 
More recently, observed solar wind plasma and 
magnetic field time series have been used to drive 
three-dimensional global MHD simulations. 
Results from these simulations have been directly 
compared to in situ spacecraft observations [e.g. 
Frank et al, 1995]. Studies of quiet 
magnetospheric intervals were considered first 
[e.g. Ashour-Abdalla et al, 1998] before studying 
the highly dynamic magnetotail during substorms 
[Raeder et al, 1998; Ashour-Abdalla et al, 1999; 
El-Alaoui, 2001]. Studies using the approach of 
using observed values as model input have 
successfully reproduced a great deal of the 
observed dynamics of the magnetosphere. These 
studies are best performed for intervals during 
which spacecraft are positioned close to regions 
where important dynamical events are occiuring. 


This paper is concerned with two substorm events 
that were first identified by Frank et al [2001]. 
Two approaches have been taken in carrying out 
comparisons between simulation results and 
spacecraft observations. In the first synthetic 
auroral emissions are derived from the simulation 
results and compared with images from auroral 
spacecraft [Fedder et al, 1995]. Because the 
aurora maps to a large region of the near-Earth 
magnetotail, this approach allows comparisons of 
multiple individual features. 

In the second approach time series from the 
simulations are compared with those measured by 
spacecraft in the magnetotail. This second 
approach has been effective in studies of different 
regions in the magnetotail. Such studies have 
addressed the dynamics of the magnetospheric 
boundary [Berchem et al, 1998] and the plasma 
sheet during substorm intervals [Ashour-Abdalla et 
al, 1999]. In these studies, the simulation results 
and magnetotail observations showed general 
agreement in their description of the spacecraft 
motion through different regions, thou^ many 
observed details were not reproduced by the 
simulations, particularly in regions with sharp 
gradients. In spite of its limitations MHD modeling 
is the only approach that can reproduce the effect 
of solar wind variations on the magnetosphere 
which then can be directly compared with satellite 
time series in specific cases. 
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2. Model 

The global MHD simulation code solves the 
resistive MHD equations as an initial value 
problem for the magnetosphere and a potential 
equation for the ionosphere \Raeder et al, 1998]. 
Since the use of a high-order hybrid scheme wifli 
increased spatial resolution minimizes the amount 
of numerical dissipation that occurs in the 
computation, we also include a resistive term in 
Ohm's law, E = -vxB + r|J;the resistivity q is a 
nonlinear function of the local current density j 
such that q = ot/'^, where a is an empirically 
determined parameter (a « 1). To avoid spurious 
dissipation we include a threshold that is also a 
function of the local normalized current density; 
this threshold is calibrated such that explicit 
resistivity is switched on only at a few grid points 
in strong current sheets. 

The solar wind magnetic field, density, 
temperature, and velocity are imposed on the 
sunward face of the simulation box; open boundary 
conditions (^Idn = 0) are assumed for all of the 
other sides of the box. The dimensions of the 
simulation box are 25 Re in the simward direction, 
300 Re along the tail and 50 Re in each transverse 
direction. With such a large simulation domain, all 
flows at the external boundaries are in the super 
magnetosonic regime, which prevents any 
information fi'om propagating back fi'om the 
boundaries and affecting the physical processes 
occurring in the simulations. A large number of 
grid points (210) are used, and equations are 
solved on a stretched Cartesian computational grid, 
which permits us to increase substantially the grid 
density in the plasma sheet region. 

The MHD model has an inner spherical boundary 
with a radius of 3 Re. This boundary is designed to 
exclude the region where the Alfven velocity 
becomes too large to use a reasonable time step. 
Since this near-Earth region is dominated by the 
terrestrial field, this procedure does not 
significantly affect the simulation. Closure of FAC 
is ensured by solving self-consistently the 
ionospheric potential equation: 

V(S.V<D) = -7||Sin0, 

where <t> denotes the ionospheric potential, Z is the 
tensor of the ionospheric conductance, is the 


FAC at the inner magnetospheric boundary 
mapped onto the polar cap, and 0 is the inclination 
of die magnetic field at the ionosphere. The 
boxmdary condition O =0 is applied at the equator. 

Our present model uses three ionization sources to 
compute the height-integrated ionospheric Hall and 
Pedersen conductivities. The first component is a 
result of solar EUV ionization and is evaluated by 
using the Moen and Brekke [1993] model. The 
second source is made up of precipitating electrons 
associated with upward FACs and is computed by 
using the relation established by Lyons et al. 
[1979]. The third contribution comes fi'om diffuse 
electron precipitation and assumes a full loss cone 
at the inner boundary. The conductances are 
calculated fi'om the precipitation parameters (mean 
energy and energy flux) by using the Robinson et 
al. [1987] empirical relation. The ionospheric 
potential is then mapped to the irmer boimdaiy 
where it is used as a boundary condition for the 
flow velocity v = -VO x B/B^. With our present 
MHD grid, the effective spatial resolution of the 
ionospheric parameters in the auroral zone is 0.5® 
in magnetic latitude and 3® in magnetic longitude. 
A detailed description of the MHD model can be 
found in Raeder et al., [1998]. 

3. Event studies 

Frank et al. [2001] conducted an investigation of 
1 12 auroral substorm events observed by the Polar 
VIS instrument in 1996. They searched for events 
where the auroral onset region could be 
magnetically mapped to a location near Geotail 
and earthward of it. Two events were found that 
satisfied their criteria; one on July 2 and one on 
September 4. In this paper we describe global 
MHD simulations of these two substorm events. 
We focus on the results of calculations of the 
power deposited into the ionosphere and 
comparisons of this to observations. 

1.1 July 1-2, 1996 event 

On July 1-2, 1996, fi'om about 1800 UT (July 1) to 
0800 UT (July 2), the WIND spacecraft was 
located 211 Re upstream of Earth in the dusk 
sector. Figure 1 shows WIND’S measurements of 
the IMF and the solar wind plasma (solid Ines) 
and the solar wind parameters at the upstream edge 
of the simulation box (dashed lines). 
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Figure 1. Solar wind parameters for July 1-2, 
1996. The first three panels show WIND 
measurements of the magnetic field components, 
then the velocity components, proton density, and 
the last panel toe thermal pressure. The dashed 
lines in toe first three panels show toe MHD 
simulation input using toe technique described in 
section 2.2 of El-Alaoui [2001]. 

In order to keep toe magnetic field divergence free 
during this event toe method described in section 
2.2 of El-Alaoui [2001] is used to model toe three 
components of toe field. This method, which is 
based on minimiun variance, does a good job in 
reproducing toe magnitude and variation of toe 
I>^. During toe whole time interval toe IMF 5y 
component was strong and directed duskward with 
occasional decreases. The IMF 5* component was 
relatively weak and drifted from toe simward 
direction to anti-sunward. The component was 
southward until 2100 UT (July 1, 1996) when it 
turned northward for about 15 minutes. It then 
turned southward again and remained southward 
until 0100 UT. The solar wind density was about 
15 to 20 cm^, which is about 3 times toe typical 
density of toe solar wind. Meanwhile, toe solar 
wind speed was slowly decreasing from about 370 
km/s to about 340 km/s. 


Figure 2 shows toe total power deposited into toe 
ionosphere from VIS observations (dashed line) 
and calculations from toe MHD simulation of toe 
power deposited into the ionosphere as a function 
of time (heavy line). The total power was 
calculated by integrating the expression below over 
a sphere at 1 /?£ . 


Power = [ 

s 


nkT 






ds 


In this expression rig and Tg are the magnetospheric 
electron density and temperature and J\\ is the 
parallel current and AO n is the field-aligned 
potential drop \Knight, 1973]. Their values were 
calculated at the inner boundary of the simulation 
at 3 Re and mapped along magnetic field lines to 
the ionosphere by using a dipole magnetic field to 
represent the internal field of the Earth. 



Figure 2. A comparison of time series fi*om VIS 
observations and results fi*om a global MHD 
simulation of the power deposited into the northern 
ionosphere. The second panel shows the power 
associated with upward Region 1 field aligned 
currents. 

The first contribution to the power comes fi*om 
diffiise electron precipitation [Kennel and 
Petschek, 1966]. The second source is 
precipitating electrons associated with upward 
FACs. We compute the mean energy eAO\\ and the 
energy flux of precipitating electrons that are 
accelerated by the parallel potential drop 
[Knight, 1973; Lyons et al., 1979]. 
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The light line in the first panel of Figure 2 shows 
the contribution of the diffuse precipitation (the 
first term in the integral). The heavy curve shows 
the total power including both diffuse and discrete 
precipitation. The dashed line shows the observed 
power deposited into the ionosphere calculated 
fi’om VIS observations. The second panel in the 
figure shows a plot of the power due to FACs (the 
second term in the integral). Only Region 1 
currents are included in this calculation, since the 
MHD model does not fully resolve the Region 2 
field-aligned currents. 

During this event the simulations suggest that the 
discrete aurora typically contributes less than 15% 
of the power, atoough after 0230 UT it seemed 
that the discrete contribution approaches the 
diffuse contribution. Both observations and 
simulations suggest that the power deposited into 
the ionosphere showed large variation versus time. 

During this whole time interval the minimum 
power deposited into the ionosphere was ~30 GW. 
The power increased to over 100 GW during the 
time of enhanced activity. Spiro et al. [1982] 
developed an empirical relationship between the 
AE index and the precipitated power that gives 
about 100 GW for major substorms. This is in 
good agreement with the simulation results. 

Another feature of this interval was small scale 
oscillations in the solar wind dynamic pressure. 
We looked for a signature of these oscillations in 
the ionospheric parameters. The parameter that 
most closely followed the oscillations was the 
power deposited into the northern ionosphere due 
to the parallel plasma flow (inJcT^V^^s). This 
parameter may reflect ULF waves in the simulated 
magnetosphere. The solid black curves in both 
panels of Figure 3 show the parallel plasma flow. 
This power is only 0.1 to 0.3% of the total power 
deposited into the ionosphere. The dashed line in 
the first panel of Figure 3 shows the solar wind 
dynamic pressure shifted by 76 minutes to account 
for the solar wind propagation from WIND 
location; in the second panel the dashed line 
represent the IMF Bj component shifted by the 
same amount. It is interesting to note that each 
small-scale peak in the solar wind dynamic 
pressure evidently corresponds to a peak in the 
calculated power (Figure 3). 



Figure 3. The solid black curves show the 
simulated power deposited into the northern 
hemisphere ionosphere due to the parallel plasma 
bulk flow. The dashed line show time shifted solar 
wind dynamic pressure (top) and IMF Bz (bottom) 
according to the scales on the right. 0200 UT is 
the time of substorm onset. 

1.2 September 4, 1996 

Solar wind plasma parameters and the IMF 
observed by the WIND spacecraft located at 
a: = 8 1 Re are shown in Figure 4. The IMF rapidly 
turned southward at 1956 UT and remained 
southward imtil ~2250 UT, when it turned 
northward. While it was southward it imderwent 
several significant changes in magnitude; at aroimd 
2130 UT the IMF B decreased in magnitude to 
nearly zero and continued to fluctuate imtil about 
2200 UT, when it increased in magnitude again. 
The IMF then changed from about -9 nT to -3 nT 
(at ~2230 UT) and to northward (at -2250 UT) in a 
step wise fashion. The IMF By component was 
dawnward throughout most of the interval except 
for a brief excursion aroimd 2050 UT. The Bx 
component showed behavior very similar to the 
IMF Bz component. The solar wind velocity was 
on average, 375 km/s throughout the interval. A 
high density of 13 cm'^ prevailed during the 
interval of southward IMF. 

Figure 5 shows, in a similar format to Figure 2, the 
total power deposited into the northern hemisphere 
ionosphere by various components of the 
precipitation. 
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Figure 4. IMF and plasma parameters from the 
WIND spacecraft for September 4, 1996. 



Figure 5. A comparison of time series from VIS 
observations (dashed line) and result from a global 
MHD simulation (solid lines). 

The contribution of the diffuse precipitation (light 
line in top panel) and the total power including 


both diffuse and discrete precipitation (black curve 
in top panel) are shown. The dashed line shows the 
observed power deposited into the ionosphere 
calculated from VIS observations. Both show large 
variations in the power deposited into the 
ionosphere. However, the simulation and the 
observations have significant differences. The 
substorm onset in the MHD simulation preceded 
the observed onset by about 15 minutes; although 
the peak energies are comparable. This may be due 
to solar wind propagation effects and limitations in 
the ionospheric model. There is also a large peak in 
the simulated power after 2230 UT that is absent in 
the observations. This discrepancy is puzzling 
because this peak is associated with a more 
strongly southward IMF at 2200 UT (at the WIND 
locations) which should increase the power. 

Small scale oscillations in the solar wind dynamic 
pressure can again be seen. The solid black curves 
in both panels of Figure 6 show the power due 
parallel plasma flow. Again this power is only 0.1 
to 0.3% of the total power deposited into the 
ionosphere. The dashed line in the first panel of 
Figure 6 shows the solar wind dynamic pressure 
shifted by 37 minutes to account for solar wind 
propagation from the WIND location; in the 
second panel the dashed line represent the IMF 
component shifted by the same amount. It is 
interesting to note that each peak in the solar wind 
dynamic pressure again seems to correspond to a 
peak in the calculated power (Figure 4). 

4. Summary 

This study used observations and MHD modeling 
together to investigate the relationship between 
solar wind variations and energy deposited into the 
ionosphere during two substorm events. During 
periods of southward IMF the MHD simulation 
showed that an increase in the power deposited 
into the ionosphere. Large scale variations in the 
simulated power usually corresponded to those in 
the observations. These variations seem to be 
correlated with the IMF Bz component. The 
simulations also revealed short time scale 
oscillations that were driven by the dynamic 
pressure fluctuations. Despite the limitation of the 
MHD model, the simulation captured important 
features of the substorms. 
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Figure 6. The black curves show the simulated 
power deposited into the northern hemisphere 
ionosphere due to the parallel plasma bulk flow. 
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